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ABSTRACT 

The  highly  stable  propagation  characteristics  of  VLF 
transmissions  make  them  an  ideal  source  of  world-wide 
navigation  data.   The  arrival  phase  of  a  VLF  signal  is 
continuously  compared  with  a  local  standard.   The  change  in 
the  phase  difference  between  the  received  signal  and  the 
local  standard  is  converted  to  relative  motion  with  respect 
to  the  VLF  transmitter. 

This  investigation  concerns  itself  with  the  development 
as  a  whole  of  the  VLF  navigation  system.   Refinements  are 
made  in  the  antenna  system  and  the  antenna  signal  processing 
In  addition,  navigation  procedures,  tables,  and  forms  are 
presented  for  system  use. 
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I.   INTRODUCTION 

Very  low  frequency  transmissions  form  the  basis  of  a 
very  accurate  long  range  navigation  system.   The  character- 
istic of  these  transmissions  that  make  them  ideal  for  such 
a  use  is  that  their  arrival  phase  at  a  point  distant  from 
the  transmitter  is  both  stable  and  predictable. 

The  VLF  wave  propagates  through  a  waveguide  formed  on 
one  side  by  the  ionosphere  and  on  the  other  side  by  the 
surface  of  the  earth.   The  phase  of  the  received  signal 
(the  carrier  frequency)  at  some  point  distant  from  the  trans- 
mitter remains  constant  as  long  as  the  propagation  path 
(waveguide)  or  the  distance  from  the  transmitter  does  not 
change.   With  the  onset  of  darkness,  the  ionosphere  height 
increases  and  modifies  the  waveguide,  causing  a  change  in  the 
observed  arrival  phase.   However,  the  observed  arrival  phase 
will  return  to  its  original  value  when  the  entire  propagation 
path  is  again  in  daylight,  and  the  waveguide  assumes  its 
initial  dimensions. 

With  the  addition  of  two  stable  oscillators,  one  at  the 
transmitter  to  serve  as  a  reference  for  the  master  oscillator, 
and  the  second  at  some  remote  receiving  site,  motion  of  the 
receiver  can  be  detected  by  observing  the  phase  of  the 
received  signal  as  compared  to  the  local  oscillator.   By 
commencing  to  track  a  transmitter  signal  in  this  manner  from 


a  known  geographic  position,  the  navigator  continuously  has 
one  line  of  position.   Proper  selection  of  a  second  trans- 
mitter will  result  in  a  second  line  of  position  crossing  the 
first,  and  an  electronic  fix. 

The  next  section  presents  the  background  of  this  system, 
and  the  basic  work  accomplished  in  its  development.   Chapter 
III  goes  on  to  present  an  explanation  of  mode  theory,  which 
describes  the  propagation  mechanism.   This  discussion  leads 
to  the  basis  of  the  system  stability,  and  hence  its  accuracy. 
Three  basic  problems  are  encountered  in  the  use  of  this  ' 
navigation  system.   These  are  examined  in  general  in  Chapter 
IV.   The  problems  are  long  path  signal  contamination,  pre- 
diction of  the  diurnal  shift,  and  the  conversion  of  phase 
information  into  geographic  positions.   Each  is  treated 
respectively  in  Chapters  V,  VI,  and  VII  with  their  solutions. 
Chapter  VIII  sets  forth  the  system,  and  treats  the  units 
from  an  individual  component  view  as  well  as  a  discussion  of 
the  system  as  a  whole.   Chapter  IX  presents  the  proposed 
procedures  for  use  of  the  system. 


II.   BACKGROUND 

The  development  of  the  long  range  VLF  navigation  system 
starts  with  Blackband  [1]  and  Stanbrough  [2],  who  both  used 
this  method  experimentally  with  good  results.   Woods  Hole 
Oceanographic  Institute  has  developed  a  table  of  geographic 
positions  [3]  which  contains  bearings  and  distances  to  various 
VLF  transmitters.   While  this  table  appears  to  be  the  initial 
effort  in  this  line,  the  geographic  positions  are  each  one 
degree  of  latitude  and  longitude  apart,  a  spacing  too  great 
for  general  navigation  use. 

Lake  [4]  conducted  a  general  investigation  of  this  VLF 
navigation  system,  covering  the  basic  propagation,  phase- 
tracking  receivers  used,  and  general  navigation  procedures. 

McKay  and  Preston  [5]  covered  in  depth  the  factors  which 
degrade   the  system  accuracy,  and  developed  the  idea  of  using 
an  antenna  with  a  cardiod-shaped  radiation  pattern  to  reject 
the  long  path  signal  which  could  contaminate  the  desired 
signal  in  some  cases. 

Roeder  [6]  developed  an  antenna  system  producing  a 
cardiod-shaped  pattern  utilizing  two  loops  and  a  whip  with 
the  associated  equipment  necessary  to  electronically  rotate 
the  resulting  pattern. 

This  paper  will  bring  together  the  efforts  of  Stan- 
brough, Lake,  McKay,  Preston,  and  Roeder  to  result  in  a  system 


for  long  range  VLF  navigation.   The  specific  results  of  this 
investigation  are  the  development  of  a  practical  shipboard 
unit  for  use  as  a  long  range  VLF  navigation  system,  an 
improved  method  for  predicting  the  time  and  amount  of  diurnal 
shift,  the  formulation  of  a  means  to  generate  navigational 
tables  of  positions,  and  the  presentation  of  specific  forms 
and  procedures  for  system  utilization. 
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III.   MODE  THEORY 

The  propagation  of  electro-magnetic  radiation  i'n  the 
VLF  range  is  best  described  by  waveguide  mode  theory.   The 
guide,  which  channels  or  contains  the  wave,  is  formed  by  the 
earth  on  the  lower  side,  and  the  ionosphere  on  the  upper 
side,  where  the  term  "ionosphere"  refers  to  the  ionospheric 
height  for  which  the  conductivity  is.  equal  to  10    mhos/m. 

As  in  conventional  waveguides ,  many  modes ,  other  than 
the  principal  mode  determined  by  the  physical  geometry, 
exist  in  the  near  region  of  sources  and  discontinuities.   The 
resulting  signal  sensed  at  any  point  in  the  field  is  equal 
to  the  vector  sum  of  the  individual  fields  due  to  the  multiple 
modes  present.   The  higher  order  modes  have  greater  attenu- 
ation factors  than  the  principal  mode,  and  therefore  their 
contributions  to  the  resultant  signal  decrease  as  the  observer 
moves  away  from  the  transmitter.   Table  I  gives  various 
values  of  attenuation  factors  for  modes  1,  2,  and  3  of  a  20 
kHz  signal  [7] . 


MODE 

ATTENUATION  FACTOR 

n  =  1 

day 
night 

1.5  dB/1000  km 

1.6  dB/1000  km 

n  =  2 

day 
night 

5.7  dB/1000  km 
2.7  dB/1000  km 

n  =  3 

day 
night 

>15   dB/1000  km 
7.5  dB/1000  km 

TABLE  I 
VLF  ATTENUATION  FACTORS 
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As  an  illustration,  examine  the  20  kHz  signal  at   1000 
km,  2000  km,  and  4000  km  from  the  source,  during  darkness. 
At  1000  km,  the  principal  mode  (n  =  1)  is  down  by  1.6  dB. 
The  second  mode  is  down  by  2.7  dB,  and  the  third  mode  is  down 
by  7.5  dB.   At  2000  km,  the  modes  are  down  respectively 
3.2  dB,  5.4  dB,  and  15  dB.   At  4000  km  the  principal  mode  is 
down  by  6.4  dB,  the  second  mode  is  down  10.8  dB,  and  the 
third  mode  is  down  30  dB.   Further  increases  in  distance 
reduce  the  signal  contributions  of  the  higher  order  modes 
and  they  may  be  neglected  from  consideration.   In  general, 
at  distances  greater  than  4000  kilometers,  the  effects  of 
higher   order  modes  (n  =  2  and  greater)  can  be  neglected. 

VLF  propagation  has  also  been  described  in  terms  of  ray 
theory.   Both  mode  theory  and  ray  theory  will  provide  the 
same  results,  however  the  mode  theory  is  better  suited  in 
this  case. 
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IV.   BASIC  PROBLEMS 

There  are  three  major  problems  to  be  solved  before  the 
system  can  be  used  for  navigation.   They  are  the  prediction 
of  the  time  of  occurrence  and  amount  of  diurnal  shift,  the 
means  of  nulling  or  canceling  the  long  path  signal  when  the 
receiver  is  in  the  multipath  zone,  and  the  problem  of  conver- 
ting phase  difference  information  into  geographic  motion 
and  position.   Each  of  these  areas  is  discussed  briefly  below, 
and  examined  in  detail  in  following  sections. 

A.   DIURNAL  SHIFT  PREDICTION 

Consider  an  example  where  the  receiving  site  is  located 
to  the  west  of  the  transmitter  and  both  places  are  in  day- 
light.  The  recorded  phase  difference  between  the  received 
signal  and  a  local  oscillator  should  have  a  constant  value 
since  there  is  no  motion  toward  or  away  from  the  transmitter, 
and  the  propagation  path  is  not  changing.   With  the  occurrence 
of  sunset  at  the  transmitter,  the  recorded  phase  difference 
at  the  receiver  starts  to  increase  as  the  height  dimension 
of  the  earth-ionosphere  waveguide  increases.   The  increasing 
height  causes  a  decrease  in  the  phase  velocity.   The  decreased 
velocity  leads  directly  to  the  change  or  increase  in  the 
observed  phase  difference.   The  phase  difference  continues 
to  increase  until  sunset  occurs  at  the  receiving  site.   This 
increase  in  the  phase  difference  is  the  diurnal  shift.   With 
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the  propagation  path  completely  in  darkness,  the  phase 
difference  is  constant  until  sunrise  occurs  at  the  trans- 
mitter.  This  time  the  phase  difference  decreases  until  it 
reaches  its  daylight  value  where  it  again  tracks  constant. 
Since  long  range  navigation  at  sea  is  a  24  hour  requirement, 
it  is  highly  desireable  that  the  amount  and  time  of  the 
diurnal  shift  be  predictable.   This  in  fact  can  be  achieved 
and  will  be  discussed  in  a  later  section. 

B.   LONG  PATH  SIGNAL 

Utilizing  the  virtual  ionosphere-earth  waveguide  and  the 
low  attenuation  factor  for  the  principal  mode,  the  signal 
radiated  from  the  antenna  has  three  prominent  characteristic 
regions,  each  being  associated  with  certain  sections  of  the 
"guide."   In  an  area  (circle  of  radius  4000  kilometers) 
around  the  transmitter,  the  multi-mode  zone  exists.   This 
region  is  not  suitable  for  phase  tracking  due  to  the  many 
modes  present.   Therefore,  this  area  is  not  useable  for 
navigating.   From  the  multi-mode  zone  outwards  from  the 
transmitter  is  the  single  mode  zone  which  is  useful  for  phase 
tracking,  and  likewise  for  navigation.   The  single  mode  zone 
extends  to  a  region  which  surrounds  the  antipode  of  the 
transmitter.   Herein  is  the  multi-path  zone.   This  region 
exists  since  the  transmitter  radiates  its  energy  omni- 
directionally.  Therefore,  in  the  area  of  the  antipode,  two 
signals  are  possible,  the  direct  or  short  path  signal,  and 
the  long  path  signal,  that  is  one  that  has  travelled  more 
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than  halfway  around  the  world.   The  presence  of  these  two 
signals  can  contribute  major  navigational  errors  if  they  are 
not  processed  with  care.   The  situation  can  go  from  a  mild 
contamination  of  the  desired  signal  to  complete  dominance 
of  the  undesired  signal.   Of  particular  concern  is  the  case 
where  minor  errors  could  be  introduced  into  the  receivers 
by  the  presence  of  this  long  path  signal,  and  these  would 
not  be  immediately  detected  by  the  navigator  since  the  system 
would  continue  to  track  in  the  proper  direction. 

The  degree  of  contamination  will  be  determined  by  the 
lengths  of  the  two  signal  paths,  and  the  signal  attenuation 
factors  for  each  path.   In  general,  seawater  paths  have 
smaller  attenuation  factors  than  land  paths ,  and  daylight 
paths  have  greater  attenuation  than  darkness  paths.   There- 
fore, a  combination  of  darkness  and  seawater  for  the  long 
path  can  provide  a  signal  of  comparable  strength  to  that  of 
the  short  path  where  it  may  have  a  significant  portion  over 
land  and  in  daylight. 

It  is  highly  desirable  to  use  the  multi-path  zone  in  that 
it  is  also  a  single  mode  zone,  and  has  very  stable  signals 
for  phase  tracking.   If  this  region  was  not  used,  the  system 
would  be  significantly  restricted.   The  use  of  an  antenna 
system  with  a  directive  pattern  which  discriminates  against 
the  undesired  signal  would  solve  this  problem. 

C.   TABLES 

In  order  to  convert  the  phase  change  (receiver  information) 
into  geographic  motion  and  position,  suitable  tables,  similar 
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to  those  existing  for  LORAN  A,  must  be  generated  to  allow 
plotting  the  ship's  position  on  charts.   These  tables  should 
include  known  geographic  positions,  and  bearings  and 
distances  (expressed  both  in  nautical  miles  and  microseconds) 
to  transmitting  stations.   The  geographic  positions  should 
be  spaced  no  greater  than  thirty  miles  apart. 

D.   MINOR  PROBLEMS 

Four  minor  problems  associated  with  maintaining  an 
accurate  plot  of  a  ship's  position  are: 

a.  Cycle  slipping 

b.  Transmitter  outage 

c.  Frequency  offset  of  the  transmitter 
master  oscillator 

d.  Sudden  ionospheric  disturbance  (SID) . 

Each  of  the  above  should  pose  no  problem  for  a  navigator 
with  a  minimum  familiarity  with  this  system.   Cycle  slipping 
refers  to  the  situation  when  the  receiver  is  tracking  a 
changing  arrival  phase  condition,  and  drops  behind  in 
indicated  phase  by  360  degrees  or  integer  multiples  thereof. 
This  would  generally  occur  only  where  the  phase  difference 
is  changing  rapidly.   This  could  only  happen  at  the  onset 
or  offset  of  the  diurnal  shift.   The  type  of  receivers  used 
in  this  investigation  allows  the  operator  to  "drive"  the 
receiver  back  onto  the  proper  cycle  if  slipping  does  occur. 
The  operator  also  has  an  option  to  choose  a  shorter  time 
constant  for  the  receiver  to  insure  that  cycle  slipping 
does  not  occur. 
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Scheduled  transmitter  outages  are  no  problem.    Normally, 
these  outages  are  scheduled,  and  the  navigator  should  have 
shifted  the  receiver  to  an  alternate  transmitter.   Tracking 
of  a  transmitter  can  commence  from  any  known  position.   By 
selecting  an  optimum  time,  the  shift  of  stations  can  be 
made  without  difficulty,  and  without  a  loss  in  position 
accuracy.   Since  a  prudent  navigator  values  a  three  lines  of 
position  fix,  risk  from  this  and  other  anomolous  outages 
can  be  minimized  by  tracking  three  transmitters  at  all  times. 

Frequency  offsets  of  transmitter  master  oscillators  are 
scheduled  sufficiently  in  advance  that  the  amount  of  the 
offset  and  date  and  time  of  its  application  will  be  promul- 
gated to  all  concerned.   The  navigator  at  sea  need  only 
make  an  appropriate  correction  for  his  calculations. 

Sudden  ionospheric  disturbances  are  caused  by  solar  flares, 
magnetic  storms,  and  high  altitude  nuclear  blasts.   The 
duration  and  magnitude  of  the  distrubance  will  depend  on  the 
duration  and  magnitude  of  the  disturbing  phenomena.   In 
general  the  disturbances  are  not  intense,  and  of  short  length. 
In  any  case,  their  effect  is  to  lower  the  ionospheric  height, 
causing  an  increase  in  the  phase  velocity  of  propagation. 
This  increase  would  make  the  transmitters  being  tracked 
appear  closer  to  the  ship.   Continuous  chart  recorders  would 
help  the  navigator  to  recognize  the  onset  of  a  distrubance. 
There  would  be  a  rapid  decrease  in  the  displayed  phase 
difference  information.   Under  this  circumstance,  a  dead 
reckoning  track  should  be  plotted  from  the  last  accurate 
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position.   Subsequent  positions  may  be  determined  by 
celestial  navigation  or  other  means.   These  positions 
compared  with  the  system  performance  will  provide  an  indica- 
tion of  when  the  ionospheric  behavior  has  returned  to  normal. 

E.   PROPAGATION  CONSIDERATIONS 

The  navigator,  in  using  this  system,  must  understand 
certain  VLF  propagation  characteristics  to  intelligently 
assess  the  validity  of  signals  under  various  circumstances. 
For  example,  it  must  be  realized  that  the  velocity  of  signal 
propagation  is  dependent  on  the  height  of  the  ionosphere 
and  also  on  the  signal  path.   The  signal  path  is  classified 
as  to  being  over  seawater,  land,  or  ice.   If  the  path  is  over 
ice,  the  ground  conductivity  is  such  that  the  wave  behaves 
as  if  it  were  propagating  over  land.   Hereafter,  reference 
will  be  made  only  to  the  land  path  due  to  the  similar 
behavior  of  propagation  over  ice.   Since  this  is  a  relative 
navigation  system,  the  percent  of  signal  path  over  land  or 
seawater  has  no  effect  on  the  system  ability  to  fix  a 
position  during  daylight  hours.   However,  as  it  will  be  seen 
in  a  following  section,  the  observed  diurnal  shift  for  the 
land  portion  of  the  path  is  not  proportional  to  that  for  the 
seawater  path.   Therefore,  the  navigator  must  be  cognizant 
of  the  amount  of  propagation  path  over  seawater,  and  over 
land  to  accurately  predict  the  diurnal  shift  corrections. 
If  a  polar  path  is  being  utilized,  he  must  also  be  knowledge- 
able of  the  extent  of  the  ice  pack  for  that  particular  time 
of  the  year. 
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As  an  example  of  the  differences  in  velocities  of  signal 
propagation,  TABLE  II  demonstrates  their  magnitudes  for  a  20 
kHz  signal  under  various  conditions.   The  four  starred 
conditions  denote  the  standard  ionospheric  conditions.   Using 


Phase 
Velocity  (m/s) 

Day/Night 

Path 

Ionosphere 
Height  (km) 

2.95 

X 

108 

Day 

Seawater 

70 

2.94 

X 

108 

Night 

Seawater 

90 

2.97 

X 

108 

Day 

Land/Ice 

70 

2.96 

X 

108 

Night 

Land/Ice 

90 

2.96 

X 

108 

Day 

Seawater 

60 

2.95 

X 

108 

Day 

Seawater 

80 

TABLE  II 
PHASE  VELOCITIES  UNDER  VARIOUS  CONDITIONS 

these  velocities  one  arrives  at  the  standard  conversion  value 
of  6.18  microseconds  per  nautical  mile  for  daytime  propagation 
over  seawater.   Night  propagation  over  seawater  has  a  conver- 
sion of  6.22  microseconds  per  nautical  mile. 

Navigation  tables  are  calculated  later  in  this  investiga- 
tion.  The  entries  are  on  the  basis  of  6.18  microseconds  per 
nautical  mile.   They  may  be  adjusted  for  the  second  conversion 
(6.22)  by  multiplying  the  distances  expressed  in  microseconds 
by  a  factor  of  1.006. 
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V.   ANTENNA  SYSTEM 

The  ability  to  discriminate  against  the  long  path  signal 
when  the  ship  is  located  in  the  multi-path  zone  is  the 
function  of  the  directive  antenna  system.   An  ideal  system 
should  have  small,  lightweight,  and  stationary  elements. 
However,  normal  multi-element  arrays  utilized  as  antennas 
in  the  VLF  range  would  be  of  such  size  that  only  a  land-based 
installation  could  support  such  configurations. 

The  use  of  a  loop  antenna,  whose  vertical  plane  is  pointed 
at  the  transmitter,  provides  a  means  of  attaining  the  above 
favorable  antenna  properties  with  bi-directional  sensitivity. 
When  this  pattern  is  coupled  with  a  properly  phased  whip 
signal,  a  cardiod-shaped  pattern  results  with  a  null  opposite 
to  the  bearing  of  maximum  signal  gain. 

The  loop  antenna  output  is  the  difference  between  the 
signals  detected  by  opposite  legs  of  the  loop.   A  signal  or 
voltage  difference  exists  due  to  the  spatial  separation 
of  the  two  antenna  legs.   If  the  antenna's  orientation  is 
changed,  the  magnitude  of  the  detected  difference  signal 
decreases  as  the  spatial  separation  (and  hence  the  phase 
difference)  with  respect  to  the  wave  front  decreases. 
Therefore,  a  null  occurs  when  the  plane  of  the  antenna  is 
perpendicular  to  the  wave  front. 

To  attain  the  cardiod  shaped  antenna  radiation  (field 
intensity)  pattern,  the  whip  antenna  signal  is  matched  in 
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magnitude  to  that  of  the  loop  antenna,  but  with  a  90  degree 
phase  lag.   The  result  of  adding  the  two  antenna  outputs  is 
a  signal  of  twice  the  strength  in  the  preferred  direction, 
and  a  null  180  degrees  from  the  maximum.   Figure  1  is  a 
composite  of  the  radiation  patterns  of  the  single  loop  antenna, 
the  whip  antenna,  and  the  resulting  cardiod  shaped  pattern. 

Since  the  cardiod-shaped  pattern  is  the  result  of  a  phase 
addition,  and  the  navigation  system  utilizes  changes  in 
phase  difference  to  detect  motion,  the  change  in  phase  across 
the  axis  of  the  pattern  is  of  major  concern.   McKay  and 
Preston  [5]  demonstrated  both  theoretically  and  experimentally 
that  a  "zero  phase  change  window"  does  exist  on  the  pattern 
axis.   This  investigation  also  confirms  the  presence  of  this 
window.   Since  this  window  is  of  finite  width,  errors  will 
not  be  introduced  in  the  system  if  the  antenna  pattern  is 
not  pointed  exactly  at  the  transmitter.   If  the  resolver 
rotor  indicates  a  direction  within  5  degrees  of  the  transmitter 
bearing,  no  appreciable  error  exists.   McKay  and  Preston  also 
have  shown  that  the  pattern  null  is  down  30  dB  for  a  sector 
of  40  degrees  width.   This  investigation  again  confirms  this 
performance.   Therefore,  any  observed  phase  changes  are 
indicative  of  the  motion  of  the  receiver  toward  or  away  from 
the  transmitter  and  not  due  to  a  phase  change  introduced  by 
the  antenna  system. 

The  use  of  a  single  loop  and  single  whip  for  phase  track- 
ing is  impractical  in  that  only  one  station  may  be  tracked 
for  each  antenna  installation,  and  also  a  sizeable  and 
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Figure  1.   Composite  Antenna  Patterns 
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massive  antenna  drive  system  would  be  required  to  continuously 
point  the  loop  antenna.   Roeder  [6]  developed  a  means  by  which 
2  stationary  loops  duplicate  the  performance  of  a  rotatable 
single  loop.   As  seen  in  Figure  2,  the  2  loops  (mounted 
concentrically  and  perpendicular  to  each  other)  drive  field 
coils  mounted  in  quadrature.   The  field  generated  in  each  coil 
is  proportional  to  the  signal  sensed  by  the  corresponding 
antenna.   The  rotor  will  sense  the  maximum  signal  when 
positioned  at  an  angle  that  corresponds  to  the  direction  of 
signal  arrival.   The  antenna  radiation  pattern  illustrated 
in  Figure  2.c  is  obtained  by  turning  the  rotor  through  36  0 
degrees.   This  antenna  system  duplicates  the  behavior  of  a 
single  loop  pointed  directly  at  the  transmitter. 

As  a  demonstration  of  the  feasibility  of  this  approach, 
Roeder  used  goniometers  in  place  of  the  field  and  sense  coils 
to  test  a  system  capable  of  tracking  two  transmitters.   A 
means  of  positioning  the  goniometer  rotors  was  required  to 
keep  the  antenna  patterns  pointed  at  the  transmitters  while 
the  ship  maneuvered.   This  was  accomplished  by  use  of  a  gyro 
compass  signal  input  to  a  power  amplifier  which  fed  a  servo 
motor  attached  to  the  rotor  shafts.   With  the  addition  of  a 
whip  antenna  and  cardiod  units,  the  proper  radiation  patterns 
were  formed. 

In  part,  this  investigation  was  directed  toward  the 
realization  of  a  more  practical  shipboard  installation.   To 
reach  this  aim,  the  six  feet  diameter  loop  antennas  used  by 
Roeder  were  replaced  by  a  submarine  VLF  crossed  loop  antenna. 
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Figure  2. a.  Basic  Schematic 


Figure  2.b.  Geographic      Figure  2.c.   Antenna 
Picture  Pattern, 

Figure  2.   Effective  Single  Loop  Antenna  System, 
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This  greatly  reduces  the  size  of  the  topside  or  mast-mounted 
unit.   An  order  of  magnitude  decrease  is  possible  since  the 
submarine  antenna  loops  are  wound  on  ferrite  cores. 

In  place  of  goniometers,  precision  resolvers  were  used 
with  a  reverse  signal  flow.   Again,  as  in  the  case  of  the 
antenna,  a  much  smaller  physical  size  is  realizable  for  an 
operational  unit,  and  as  a  result  allows  the  direct  use  of  a 
gyro  repeater  signal  for  driving  or  positioning  the  rotors 
without  the  need  of  a  power  amplifier. 

The  gyro  compass  repeater  input  can  be  either  120  VAC, 
60Hz  or  120  VAC,  400  Hz.   The  synchro  receiver  drives  a  local 
compass  card.   Two  other  compass  cards  are  geared  to  the 
driven   unit.   A  means  of  securing  the  rotors  of  the  resolvers 
to  the  compass  cards  is  provided  through  a  bearing  pointer. 
The  bearing  of  the  station  to  be  tracked  is  indicated  by  the 
pointer  when  the  unit  is  in  operation. 

Audio  amplifiers  are  used  to  boost  the  VLF  signal  level 
for  the  two  loops'  inputs  into  the  resolvers.   The  two 
amplifiers  feed  four  windings.   Each  amplifier  feed  corre- 
sponding windings,  one  in  each  resolver.   Separate  gain 
controls  are  present  for  each  winding  in  each  resolver.   More 
resolvers  may  be  fed  to  allow  a  capability  to  simultaneously 
track  more  transmitters.   Caution  should  be  exercised  to 
insure  that  sufficient  signal  energy  is  available  at  the 
output  of  the  amplifiers  for  this  purpose. 

In  order  to  combine  the  loop  and  whip  signals  properly, 
commercial  cardiod  units  are  used  in  this  system.   These 
units  are  manufactured  by  TRACOR,  Inc.  of  Austin,  Texas. 
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They  provide  a  means  to  vary  independently  the  magnitude  and 
phase  of  the  whip  signal.   The  loop  signal  may  be  reversed 
in  phase.   The  output  of  the  two  cardiod  units  furnish  the 
input  or  antenna  signals  for  the  VLE1  tracking  receivers. 

Figures  3,  4,  and  5  represent  the  results  of  tests 
performed  on  the  antenna  system  using  VLF  transmitters  NAA, 
NBA,  and  NPM  as  radiating  sources.   Each  of  the  figures 
contain  two  plots.   A  circular  plot  of  the  field  intensity 
demonstrates  in  each  case  that  the  axis  of  the  pattern  points 
in  the  direction  of  the  transmitter.   Also,  the  cardiod-shaped 
patterns  are  well  defined,  and  clearly  exhibit  the  null  and 
severe  attenuation  in  the  reverse  direction.   The  second  plot 
is  the  measure  of  the  phase  distortion  error.   These  plots 
are  centered  on  the  bearing  of  the  transmitter  and  span  a 
sector  of  20  degrees.   These  measurements  confirm  the  existance 
of  the  "zero  phase  change  window"  on  the  pattern  axis,  and 
the  proper  formation  of  the  antenna  pattern. 

The  tests  were  conducted  in  the  following  manner.   Two 

receivers  were  used  to  track  a  transmitter.   One  receiver  used 

the  antenna  system  as  a  signal  source,  and  the  second  receiver 

used  a  whip  antenna  only.   This  second  receiver  served  as  a 

reference  for  comparison  with  the  performance  of  the  antenna 

system  channel.   The  field  intensity  was  measured  every  5 

degrees  of  rotation  of  the  resolver  rotor,  and  plotted  as 

shown  in  Figures  3. a,  4. a,  and  5. a.   For  the  phase  distortion 

error  measurements,  the  phase  difference  data  of  both  receivers 

was  compared  as  the  resolver  rotor  position  was  changed.   This 

information  is  presented  for  a  20  degree  window  in  Figures  3.b, 

4.b,  and  5.b. 
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Figure  3. a.   Intensity  Pattern. 
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Figure  3.   Antenna  Field  Intensity  Pattern  and 
Phase  Distortion  Error  Source  -  NBA, 
Bearing  115. 
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Figure  4.   Antenna  Field  Intensity  Pattern 
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Figure  5.   Antenna  Field  Intensity  Pattern  and 
Phase  Distortion  Error  Source  -  NPM, 
Bearing  255. 
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VI.   DIURNAL  SHIFT 

The  observed  relative  phase  difference  between  the  stable 
local  oscillator  and  the  received  VLF  signal  (hereafter 
referred  to  as  the  "phase  difference")  at  a  remote  receiver 
increases  as  the  ionospheric  height  increases  with  the 
occurrence  of  sunset.   The  phase  difference  settles  at  a  new 
constant  value  when  the  complete  path  from  the  transmitter 
to  the  receiver  is  in  darkness.   The  reverse  process  of 
sunset  is  observed  at  sunrise.   This  change  between  the  two 
constant  values  of  the  phase  difference  is  called  the  diurnal 
shift.   In  order  that  the  navigator  may  determine  his  position 
from  sunset  to  sunrise,  the  diurnal  shift  must  be  predictable 
in  the  time  of  occurrence/  and  the  amount  of  additional  phase 
difference. 

McKay  and  Preston  [5]  approached  this  problem  by  calcu- 
lating various  relative  phase  velocities  for  different 
conditions.   A  more  direct  solution  to  this  problem  is  a 
combination  of  the  work  of  Wait  [8]  and  Brady  and  Crombie  [9]. 
Wait's  work,  which  is  described  in  the  next  paragraph,  gives 
us  the  amount  of  diurnal  shift,  and  Brady  and  Crombie  provide 
the  means  to  calculate  the  time  of  sunset  and  sunrise  at  the 
height  of  the  ionosphere.   The  work  of  Brady  and  Crombie  will 
be  treated  more  fully  in  a  following  section  on  navigation 
procedures . 

This  solution  assumes  that  the  daily  phase  difference 
variation  for  a  stationary  receiving  site  can  be  accurately 
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represented  by  a  trapezoid  as  shown  in  Figure  6 .   For  this 
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Figure  6.   Diurnal  Phase  Shift. 

example,  consider  the  case  where  the  receiving  site  is  to  the 
west  of  the  transmitter.   During  the  period  represented  by 
the  interval  A  to  B,  there  is  a  constant  phase  difference 
observed.   At  the  time  B,  sunset  occurs  at  the  transmitter, 
and  the  height  of  the  ionosphere  starts  to  increase  causing 
an  increase  in  the  phase  difference  due  to  the  decreasing 
phase  velocity.   The  phase  difference  continues  to  increase 
until  sunset  occurs  at  the  receiver.   At  this  time,  the 
waveguide  path  has  again  stabilized,  and  a  constant  phase 
difference  is  observed  until  sunrise  occurs  at  the  transmitter 
at  time  D.   Then  the  phase  difference  decreases  until  time  E 
when  sunrise  occurs  at  the  receiver.   The  phase  difference 
is  again  at  the  same  value  as  for  the  preceding  day.   If  the 
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receiver  was  to  the  east  of  the  transmitter,  point  B  would 
represent  the  time  of  sunset  at  the  receiver  and  C,  the  time 
of  sunset  at  the  transmitter.   Likewise,  D  would  be  the  time 
of  sunrise  at  the  receiver,  and  E  would  be  the  time  of  sunrise 
at  the  transmitter.   However,  the  measured  value  of  At  would 
be  the  same  in  both  cases. 

Wait  developed  the  following  equation 

At  =  ~  (d/v   -  c/v,) 
.  3    '  n    '    d 

where  At  is  the  total  diurnal  phase  shift  in  microseconds, 
D  is  the  distance  in  kilometers  between  the  transmitter  and 
receiver,  c  is  the  velocity  of  light,  v   is  the  phase  velocity 
of  the  signal  at  night,  and  v,  is  the  phase  velocity  during 
the  day,  the  last  three  quantities  having  units  of  meters  per 
second. 

The  foregoing  presupposed  a  stationary  receiving  site. 
Now,  with  the  receiver  located  on  a  ship  that  is  moving  away 
from  the  transmitter,  Figure  6  must  be  modified  by  adding 
point  by  point  a  second  plot  which  is  a  line  with  a  slope 
equal  to  N  x  6.18  microseconds  per  hour  where  N  is  equal  to 
the  speed  of  the  ship  in  knots  in  the  direction  of  the 
transmitter.   (Note,  the  slope  is  positive  for  motion  away 
from  the  transmitter,  and  negative  for  motion  toward  the 
transmitter.)   At  the  same  time,  the  night  or  darkness  portion 
of  the  trapezoid  pattern  will  be  modified  since  At  varies 
linearly  with  D.   Therefore,  for  a  more  complete  picture  of 
the  diurnal  shift  as  it  applies  basically  to  the  VLF  navigation 
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system,  Figure  7  is  presented  with  the  various  portions  of 
the  composite  process. 
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Figure  7.   Diurnal  Shift  with  Ship  Motion. 


Examining  Wait's  equation,  the  distance  to  the  transmitter 
must  be  known,  and  ratios  of  the  velocity  of  light  to  the 
phase  velocity  for  both  day  and  night  must  be  determined. 
The  great  circle  distance  can  be  calculated  using  methods 
such  as  Bowditch  [10].   In  order  to  find  the  velocity  ratios, 
an  equation  from  the  Omega  System  handbook  [11]  is  used: 
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v  h. 

-£=  1  -  0.36  —  + 
c  a 


(27m  -  V      -<P.)  £ 

g     1   4tt/2  f  h. 


=   ionospheric  height  for  which  the  conductivity 
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Where 

v  =  signal  phase  velocity 

c  =  free  space  velocity  of  light 

0.36  =  constant  chosen  to  give  a  good  fit 

h 

i 

is  equal  to  10  w  mhos/m 
=   radius  of  the  earth  =  6378.155  kilometers 
=   phase  shift  on  reflection  at  the  ground 
=   phase  shift  on  reflection  at  the  ionosphere 
=   mode  number  (n  -  1  for  the  principal  mode) 
f    =   frequency  in  Hz 

Initially,  an  attempt  was  made  to  compare  the  results 
of  the  observed  diurnal  shift  with  the  calculated  value  of 
At  using  the  above  equations.   Since  the  distance  was  fixed, 
it  was  possible  to  arrive  at  accurate  values  for  (c/v 
c/vj .   The  measured  values  did  not  agree  with  those 
calculated. 

A  second  attempt  to  predict  the  diurnal  shift  began  with 
measuring  At  for  stations  where  the  propagation  path  was 
entirely  over  land  (stations  NAA,  NSS ,  and  NBA).   This  resulted 
in  three  accurate  values  for  (c/v   -  c/v,) .    The  heights  used 
for  the  ionosphere  were  70  kilometers  for  day,  and  90  kilo- 
meters for  night.   The  constant,  0.36,  was  replaced  by  "C" , 
and  the  equation  solved  for  "C".   This  resulted  in  three 
values  for  "C"  as  shown  in  Figure  8.   As  an  initial  interpre- 
tation, "C"  was  represented  by  a  linear  function  of  frequency. 
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Therefore,  the  constant  for  a  good  fit  has  been  replaced  by 
a  linear  function  of  frequency.   Using  the  data  from  Figure 
8/  a  curve  showing  the  value  of  C   versus  frequency  is 
plotted  as  Figure  9  for  the  land  propagation  route,  where  C 
is  defined  as  (0.53961/0.3) (c/v   -  c/v.) .   This  substitution 
allows  one  to  use  nautical  miles  directly  for  the  computation 
of  At. 

The  same  procedure  was  followed  for  sea  propagation 
routes.   These  results  are  shown  in  Figures  10  and  11.   A 
similar  curve  is  shown,  Figure  12,  which  was  derived  from 
data  presented  by  Spies  and  Wait  [7] .   Figure  13  is  an  overlay 
of  Figures  9,  11,  and  12  to  compare  the  general  shapes  and 
slopes  of  the  three  curves. 

Figures  9  and  11  along  with  a  modified  version  of  Wait's 

equation,  would  be  used  by  the  navigator  in  his  prediction 

of 

At  =  D  C 
v 

the  amount  of  the  diurnal  shift.   In  the  case  where  the  entire 
propagation  path  is  over  seawater,  C   is  read  from  Figure  11 
for  the  proper  frequency.   Then,  At  would  be  equal  to  the 
product  of  C   and  D.   If  the  propagation  path  for  the  signal 
is  over  both  land  and  seawater.  At  is  then  the  sum  of  the 
diurnal  shift  for  the  land  portion  of  the  path,  and  the 
seawater  portion  of  the  path.   For  example,  examine  the 
situation  where  a  ship  is  6000  miles  from  the  transmitter 
(f  =  22.0  kHz).   The  seawater  portion  of  the  path  is  4000 
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nautical  miles.   The  diurnal  shift  due  to  the  land  portion 
is  equal  to  2000  x  0.0133  or  26.6   microseconds.   The  diurnal 
shift  due  to  the  seawater  portion  is  equal  to  4000  x  0.0196 
or  78.4  microseconds.   The  total  diurnal  shift  for  this  path 
is  105.0  microseconds. 
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VII.   TABLES  OF  POSITIONS 

In  order  to  use  this  system  of  navigation,  tables  of 
positions  must  be  generated  that  will  allow  the  navigator  to 
convert  the  measured  changes  in  signal  phase  to  geographic 
coordinates . 

Since  this  system  is  based  on  relative  changes  in  signal 
phase,  some  standard  procedure  must  be  adopted  for  the 
navigational  tables.   This  investigation  used  the  great 
circle  distance  in  nautical  miles  from  a  known  position  to 
the  transmitter.   (This  distance  is  also  listed  in  micro- 
seconds to  facilitate  the  plotting  of  new  positions.   The 
generated  tables  also  list  the  bearing  of  the  transmitter 
and  a  second  bearing,  90  degrees  from  the  first  to  aid  the 
navigator  in  laying  down  lines  of  position  on  charts  and 
plotting  sheets.) 

The  known  points  are  spaced  every  fifteen  minutes  of 
latitude  and  longitude  to  reduce  the  need  for  interpolation 
over  large  distances. 

Contained  in  the  Computer  Output  section  is  a  sample  page 
from  the  tables.   Across  the  top  of  the  page  is  the  name  of 
the  transmitting  station,  its  frequency,  and  geographic 
location.   Latitude  is  given  in  the  left  most  and  right  most 
columns,  and  longitude  is  denoted  in  the  first  row.   For  any 
given  position,  there  are  four  entires.   The  first  line  is 
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the  great  circle  distance  in  nautical  miles  to  the  transmitter, 
The  second  line  contains  the  same  distance  expressed  in  micro- 
seconds.  The  third  line  contains  two  entries,  the  first  is 
the  bearing  to  the  station  along  the  great  circle  path,  and 
the  second  is  a  bearing  to  be  used  to  lay  down  a  very  good 
approximation  to  the  circular  line  of  position.   The  navigator 
will  be  concerned  with  the  bearing  of  the  transmitter  along 
the  great  circle  path,  in  order  to  have  the  "zero  phase  change 
window"  of  the  cardiod  shaped  radiation  pattern  pointed  in 
the  proper  direction.   Also,  when  the  ship  is  located  in  the 
multi-path  zone,  knowledge  of  this  bearing  will  allow  the 
radiation  pattern  null  to  be  pointed  at  the  long  path  or 
undesired  path  signal. 

The  method  of  calculation  of  the  tables  is  a  FORTRAN  IV 
program,  contained  in  the  Computer  Program  section  which  was 
run  on  an  XDS  9  300  computer  located  at  the  Naval  Postgraduate 
School.   The  mathematical  model  was  derived  from  Ref.  [10]. 
The  program  computes  two  quantities,  the  distance  and  true 
bearing  from  a  known  geographic  position  to  a  transmitter. 
The  calculations  are  based  on  the  following  two  equations, 
hav  D  =  hav  DLo  cos  L,  cos  L?  +  hav  1 
hav  C  =  sec  L,  esc  D  (hav  coL~  -  hav (D~coL, ) ) 

The  function  hav  (haversine)  of  an  angle  is  equal  to  the 

21 

sine  squared  of  half  the  angle  (hav  A  =  sin  2A) .  The  quan- 
tities in  the  equations  are  as  indicated  below.  Figure  14 
shows  the  basic  navigational  triangle  as  it  appears  in  this 
case. 
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Station 


Figure  14.   Solid  Geometry 

Navigational  Triangle 
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D    -   distance  expressed  in  degrees  between  the  known 
position  and  the  transmitter 

DLo  -   difference  in  longitude  measured  from  the  known 
position  to  the  transmitter,  and  identified  as 
either  east  or  west  (This  quantity  will  be  less 
than  180  degrees.). 

L,   -   latitude  of  the  known  position 

L«   -   latitude  of  the  transmitter 

1    -   difference  in  latitude  (1  is  measured  from  the 
known  position  latitude  to  the  transmitter 
latitude  and  identified  as  either  north  or  south.). 

C       an  angle  indicating  the  direction  of  the  trans- 
mitter from  the  known  position  labeled  north  or 
south  to  agree  with  L,  and  east  or  west  to  agree 
with  DLo 

CoL^-   colatitude  of  the  transmitter  (CoL~  =  90  -  L~  are 
of  the  same  name,  and  CoL~  =  90  +  L-  if  L,  and 
Lp  are  of  contrary  names) 
D~coL-.  -   the  numerical  difference  between  D  and  coL, 

The  computer  program  processes  the  angular  position  data 
in  two  formats,  degrees  and  radians.   The  radian  format  is 
necessary  for  use  in  the  computer  angular  subroutines.   The 
distance  is  first  calculated  and  then  used  as  an  input  for 
arriving  at  the  bearing  to  the  transmitter.   The  program 
includes  the  necessary  logic  for  determining  the  correct 
procedures  and  labels  in  the  computation  of  DLo,  1,  CoL2 ,  plus 
the  conversion  of  C  into  a  true  bearing. 
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For  the  purposes  of  this  investigation,  tables  were  . 
generated  for  an  area  of  the  Pacific  Ocean  adjacent  to  the 
coast  of  California,  centered  around  Monterey.  (The  coordinates 
of  this  area  are  from  34N  to  40N,  and  from  the  coast  out  to 
122W.)   These  tables  are  in  the  same  format  as  that   in  the 
Computer  Output  section,  and  are  in  the  custody  of  the  Dean 
of  Research  Administration  at  the  Naval  Postgraduate  School, 
Monterey,  California.   Similar  tables  may  be  produced  for 
any  other  area  of  interest  by  using  the  proper  geographic 
locations  as  input  data  for  the  computer  program. 
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VIII.   VLF  NAVIGATION  SYSTEM 

Herein  is  a  listing  of  the  navigation  system  components 
and  a  discussion  of  their  properties  and  characteristics 
as  they  apply  to  the  system.   The  units  are  referred  to  by 
their  manufacturers  and  model  numbers.   These  particular 
components  were  used  in  this  investigation  since  they  were 
available  at  the  Naval  Postgraduate  School.   Similar  units 
from  other  manufacturers  meet  the  requirements  for  this 
system  and  will  provide  the  needed  performance. 

A.   COMPONENTS 

The  basic  shipboard  system  consists  of  the  following  items 

One  stable  oscillator 

Three  VLF  phase  tracking  receivers 

One  VLF  whip  (sense)  antenna 

One  VLF  crossed  loop  antenna 

One  antenna  control  unit 

1.   Stable  Oscillator 

The  very  heart  of  the  shipboard  system  is  the  stable 
oscillator.  The  stability  of  this  oscillator  will  provide 
the  mark  and  measure  of  the  system  accuracy.  The  transmitters 
listed  in  Table  III  all  have  a  frequency  stability  of  1  part 
in  10  (The  meaning  of  "1  part  in  10  "  is  that  for  a  given 
oscillator  or  standard,  the  output  frequency  will  deviate  at 
the  most  1  hertz  in  10   hertz  from  the  nominal  value)  [12] . 
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As  an  illustration  of  system  accuracy,  let  us  assume  the 
shipboard  oscillator  has  a  stability  also  of  1  part  in  10 
Now  given  the  worst  instance  where  these  errors  are  additive, 
2  parts  in  10   ,  examine  the  system  accuracy  at  this  point. 
Since  there  are  0.86  4  x  10    microseconds  in  one  day,  the 
above  error  of  2  parts  corresponds  to  2.32  microseconds. 
Using  the  conversion  of  6.18  microseconds  per  nautical  mile, 
and  2000  yards  per  nautical  mile,  the  maximum  error  of  2  parts 
corresponds  to  750  yards.   If  the  shipboard  oscillator  has  a 
stability  of  2  parts  in  10   ,  the  corresponding  maximum 
error  increases  to  1125  yards. 

Present  day  frequency  standards  fall  into  two  general 
categories,  either  a  high  precision  temperature  controlled 
crystal,  or  a  device  utilizing  a  natural  atomic  resonance. 
There  are  two  stable  oscillators  located  at  the  Naval  Post- 
graduate School,  one  each  of  the  types  mentioned  above,  both 
of  which  are  suitable  for  system  use. 

The  high  quality  quartz  crystal  oscillator  is  a  Model 
2.5,  manufactured  by  the  Sulzer  Laboratories  of  Rockville, 
Maryland.   This  unit  requires  a  22-32  VDC,  0.33  ampere  input, 
with  outputs  (1  volt  to  50  ohms)  at  frequencies  of  2.5  mHz , 
1  mHz,  and  100  kHz.   The  crystal  unit  operates  at  2.5  mHz  and 
is  kept  at  a  constant  temperature  by  means  of  a  two  stage 
oven  [13].   The  crystal  will  drift  with  age,  and  monitoring 
of  its  output  is  necessary  to  determine  the  amount  of  drift. 

Prior  to  use  of  this  type  of  crystal  oscillator,  it  is 
necessary  that  the  unit  be  in  operation  for  periods  from  3 
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to  6  months  for  settling.   Loss  of  power  to  the  unit  will 
require  the  long  stabilization  procedure.   Therefore,  this 
type  of  oscillator  is  generally  powered  from  wet  cell  batter- 
ies, which  are  constantly  on  charge  from  a  normal  power  supply. 
On  loss  of  power,  the  oscillator  is  then  completely  carried 
by  the  batteries . 

Atomic  resonance  type  standards  utilize  the  emitted 
energy  of  a  given  natural  decay  transition  as  a  basis  of 
frequency  definition.   The  Model  V-4700  Rubidium  vapor  fre- 
quency standard  manufactured  by  Varian  Associates  of  Palo  Alto, 
California,  uses  a  hyperfine  transition  of  Rubidium  87  to 
stabilize  a  5  mHz  crystal  oscillator  [14].   This  unit  has 
outputs  of  5  mHz,  1  mHz,  and  100  kHz.   It  has  a  long  term 
stability  of  5  parts  in  10   ,  and  a  short  term  stability  of 
1  part  in  10   .   The  Varian  model  uses  a  similar  power  supply 
as  the  Sulzer  crystal  unit.   The  atomic  resonance  type 
standards  require  only  a  few  days  for  the  stabilization 
process  as  compared  to  the  months  required  for  crystal  units. 

The  Sulzer  crystal  has  a  short  term  stability  of  8 
parts  in  10   ,  which  would  correspond  to  a  maximum  error  of 
3375  yards.   This  amount  of  error  is  greater  than  an  amount 
to  be  tolerated  for  the  system.   On  this  basis,  a  stability 
of  4  parts  in  10    would  provide  accuracy  of  less  than  one 
nautical  mile. 

In  the  selection  of  a  frequency  standard,  one  has  the 
option  of  either  type  of  the  two  discussed  above  or  perhaps 
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a  Cesium  standard.   Two  performance  characteristics  must  be 
evaluated,  stability  (or  inversely  dispersion)  and  drift. 
The  Cesium  standard  has  no  drift  but  a  moderate  dispersion. 
The  dispersion  problem  could  be  remedied  by  suitable  averaging 
(integration)  techniques,  but  the  prime  drawback  for  use  in 
this  navigation  system  is  its  high  cost.   Going  to  a  Rubidium 
standard,  the  dispersion  factor  is  improved  (reduced),  but 
drift  now  enters  the  picture.   This  adds  the  requirement  for 
continual  correction  to  maintain  an  accurate  standard.   Use 
of  a  highly  precise  crystal  has  the  same  problems  as  the 
Rubidium  standard,  but  with  greater  drift. 

Therefore,  with  both  of  the  practical  standards,  drift 
must  be  determined  by  comparison  with  a  refined  standard. 
Depending  on  the  amount  of  drift  and  the  smallest  correction 
that  may  be  made  to  the  particular  frequency  standard,  a 
regular  correction  procedure  must  be  conducted  to  maintain 
an  accurate  navigation  system. 
2 .   VLF  Tracking  Receivers 

TRACOR  Series  599  VLF  Tracking  Receivers,  manufactured 
by  TRACOR,  Incorporated  of  Austin,  Texas,  were  used  in  this 
investigation.   These  receivers  were  designed  for  time  and 
frequency  measurements  with  the  VLF  transmitting  stations 
listed  in  Table  III. 

Figure  15  is  a  simplified  block  diagram  that  depicts 
the  operation  of  phase  tracking  VLF  receivers.   As  illustrated, 
the  output  of  the  phase  comparator  is  the  input  of  the 
electronic  phase  shifter.   The  circuit  is  such  that  the  phase 
error  is  driven  to  zero,  and  the  system  tracks  an  internal  null 
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This  receiver  is  powered  either  from  120  VAC (48-400 
Hz)  or  from  a  suitable  DC  source.   Useful  outputs  of  this 
receiver  are  as  follows : 

a)  time  difference  register  -  a  digital  counter  that 
indicates  the  relative  phase  difference  in  microseconds 
between  the  received  VLF  carrier  and  the  local  frequency 
standard. 

b)  audio  -  1  kHz  intermediate  frequency. 

c)  phase  difference  (100  microsecond  scale)  for  display 
on  an  external  chart  recorder. 

d)  relative  carrier  level  from  a  front  panel  meter  and 
also  from  a  recorder  output. 

The  receiver  requires  a  0.01  microvolt  signal  at  its  terminals 
to  enable  normal  phase  tracking  [15]. 
3 .   Antennas 

The  sensing  apparatus  for  this  navigation  system 
consists  of  three  separate  antennas.   Two  loop  antennas  are 
concentrically  mounted  at  right  angles  in  a  submarine  VLF 
antenna.   This  antenna  is  designated  as  an  AT-317C/BRR  (fre- 
quency range  of  14.6  -  38.0  kHz).   For  a  shipboard  installa- 
tion, one  loop  is  aligned  with  the  ship's  centerline,  and  the 
other  loop  is  at  a  right  angle  or  running  port  and  starboard. 
Each  of  these  loops  has  a  separate  but  identical  cable  run 
to  the  Antenna  Control  Unit.   A  VLF  whip  antenna  (TRACOR  Model 
No.  599  800)  is  also  used  to  provide  a  sense  input  to  obtain 
the  desired  antenna  radiation  pattern.   The  only  requirement 
placed  on  the  location  of  the  whip  antenna  relative  to  the 
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loop  antennas,  is  that  the  distance  should  be  small  compared 
to  a  wavelength.   This  requirement  provides  no  restriction 
when  operating  at  VLF  frequencies. 
4 .   Antenna  Control  Unit 

The  Antenna  Control  Unit  consists  of  the  Antenna 
Control  Console,  the  two  cardiod  units,  plus  the  associated 
auxiliary   components.   Its  purpose  is  to  combine  the  signals 
from  the  three  antennas  into  an  effective  single  antenna 
signal.   The  resulting  antenna  pattern  is  to  be  rotatable  to 
allow  the  ship  freedom  of  motion  to  maneuver,  and  also  to 
allow  the  navigator  the  freedom  of  choice  to  change  the 
station  being  tracked.   Therefore,  a  gyro  compass  repeater 
signal  is  utilized  to  provide  the  direction  sense,  and  an 
antenna  pattern  control  capability.   An  earlier  section 
discussed  the  actual  operation  and  performance  of  the 
components  and  the  unit  as  whole. 

B.   SYSTEM 

Figure  16  depicts  the  components  of  the  navigation  system 
and  their  interconnection.   Two  rotatable  antenna  radiation 
patterns  (formed  by  means  of  the  three  antennas,  the  resolvers , 
and  the  cardiod  units),  sense  VLF  signals.   These  detected 
signals  are  compared  in  the  phase-traking  receivers  with  the 
local  standard  through  the  use  of  an  internal  frequency 
synthesizer  and  a  phase-lock  loop  technique.   The  phase 
difference  (a  relative  quantity)  is  read  out  on  a  digital 
indicator  on  the  face  of  the  receiver,  and/or  displayed  on  a 
continuous  graph  recorder  as  a  function  of  time.   The  output 
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format  of  the  phase  difference  is  in  microseconds  of  change, 
which  may  be  converted  to  distance  traveled  by  the  conversion 
unit  of  6.18  microseconds  per  nautical  mile,  or  through  the 
use  of  tables. 

While  Figure  16  shows  RCVR  #3  being  fed  from  the  whip 
antenna,  a  third  resolver  channel  can  be  used  to  increase  the 
system  flexability.   The  whip  signal  is  useful  only  in  the 
single  mode  for  a  transmitter. 

The  antenna  radiation  pattern  must  be  capable  of  tracking 
the  transmitter  to  insure  that  the  received  VLF  signal  enters 
through  the  zero  phase  change  window  on  its  axis.   To  achieve 
this  capability,  it  is  best  to  consider  the  components  of  the 
antenna  patterns.   The  cardiod-shaped  pattern  is  composed  of 
the  "figure  eight"  pattern  from  the  "single  effective  loop" 
and  the  whip  pattern.   The  cardiod  orientation  is  the  same 
as  that  of  the  loop  antenna  pattern,  and  this  is  controlled 
by  the  position  of  the  resolver  rotor.   Therefore,  the  means 
of  rotating  the  resultant  cardiod-shaped  pattern  is  to  rotate 
the  resolver  rotor. 

As  course  changes  are  made,  the  ship  "rotates"  under  the 
rotor  of  the  gyro  compass,  which  maintains  its  orientation. 
Likewise,  the  true  bearing  to  a  transmitter  will  not  change 
with  the  ship's  course  change  (the  relative  bearing  measured 
from  the  ship's  head  or  bow  will  change  an  amount  equal  to  the 
course  change).   Therefore,  by  controlling  the  position  of  the 
resolver  rotor  with  a  gyro  repeater  input,  the  antenna  pattern 
will  remain  pointed  at  the  transmitter  regardless  of  the 
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ship's  maneuvers.   Small  corrections  will  be  made  as  the  ship 
traverses  large  distances,  and  the  bearings  to  each  trans- 
metter  changes.   Tables  of  known  positions,  discussed  in  the 
next  section,  will  contain  these  bearings.   The  navigator  need 
only  make  these  changes  to  the  resolver  indicators. 
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IX.   NAVIGATIONAL  FORMS  AND  PROCEDURES 

A,   FORMS 

Four  new  navigational  forms  are  proposed  for  use  with  the 
VLF  navigation  system.   Their  purposes  and  use  are  explained 
in  the  following  sections. 

Reference  Position  Form 

This  form,  Figure  17,  records  a  position  that  is  fixed 
by  outside  means.   The  form  has  a  space  to  record  the  oscilla- 
tor drift  in  microseconds  per  day.   The  entry  would  be  that 
determined  by  the  last  calibration  check.   Space  is  provided 
for  a  description  of  the  method  of  fixing  the  ship's  position, 
also.   Typically,  at  the  start  of  a  voyage,  an  appropriate 
location  can  be  used  as  a  starting  reference.   At  the  time 
of  the  navigational  fix,  the  receiver  microsecond  output  is 
recorded.   Space  is  provided  for  the  recording  of  the  geo- 
graphic line  of  position  in  microseconds  from  the  table  of 
positions.   A  is  the  numerical  difference  in  the  receiver 
output  and  the  table  of  position  value,  such  that  A  may  be 
added  to  the  receiver  reading  (also  plus  corrections)  to 
obtain  the  proper  chart  value  for  plotting  purposes. 

Hourly  Navigation  Record 

This  form,  Figure  18,  provides  the  main  log  for  hourly 
recording  of  the  ship's  location.   The  receiver  output  is 
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DATE 
TIME 

POSITION 

RCVR  #1 
STATION 

RCVR  #2 
STATION 

RCVR  #3 
STATION 

RCVR  RDG 
TABLE 
A 
OFFSET 

RCVR  RDG 
TABLE 
A 
OFFSET 

RCVR  RDG 
TABLE 
A 
OFFSET 

OSCILLATOR  DRIFT 

DESCRIPTION  OF  FIX 


FIGURE  17.   Reference  Position  Form, 
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entered  in  the  column  labelled  "RCVR".   The  total  correction 
to  the  receiver  output  is  entered  in  the  next  column  labelled 
"CORRT" .   This  should  include  A,  oscillator  drift,  and  At 
if  appropriate  for  that  particular  time  of  the  day.   The  hourly 
position  may  be  entered  in  the  two  extreme  right  hand  columns. 
Sufficient  space  is  provided  at  the  bottom  of  the  form  to 
allow  the  navigator  to  record  the  predicted  values  for  At 
for  night  navigation. 
Sunset/Sunrise  Calculation  Forms 

These  two  forms,  Figures  19  and  20,  are  used  to  calculate 
the  time  of  sunset/sunrise  at  ionospheric  heights.   Both  forms 
are  self  explanatory,  and  are  based  on  the  work  of  Brady  and 
Crombie  [9],   Appendix  A  is  a  collection  of  useful  data  for 
use  with  these  two  forms. 

B.   PROCEDURES 

Navigation  procedures  can  best  be  described  in  a  narrative 
type  discussion  of  a  navigator's  use  of  this  sytem. 

Prior  to  the  actual  use  of  the  navigation  system,  the 
navigator   selects  two  or  more  transmitting  stations  for 
tracking.   Among  the  things  to  be  considered  in  his  selection 
are  bearings  to  the  stations  in  the  area  where  the  ship  will 
be  steaming,  scheduled  transmitter  maintenance  periods,  trans- 
mitter power,  and  distance  to  the  station.   Bearings  and 
distance  to  the  transmitting  stations  will  be  listed  in  the 
tables  of  position.   Stations  should  be  selected  so  that  their, 
bearings  will  cross  at  least  30  degree  angles  (preferrable  as 
close  to  90  degrees  as  possible)  to  insure  maximal  accuracy. 
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SUNSET 
DATE 


L_ 

X_ 
6 


GHA(1200) 

SIN  L 

COS  L 

SIN  6 

COS  6 


COS (LHA)=-. 1103  -  A  x  B 


C  x  D 


LHA  =  (LHA) 


15 

GHA  <  180 

E  =  -GHA 
15 


eD  =  A_ 

U    15 


or   GHA  >  180 

E  =  360-GHA 
15 


TS 


_A 

_C 

_B 

D 


+ 

~* 

12.0000 

xxxxxxx 

E 

e    w 

6D 

E 

LHA 

XXXXXXX 

L_ 
X_ 

6 


GHA(1200) 


SIN  L 

A 

COS  L 

C 

SIN  6 

B 

COS  6 

D 

COS (LHA)=- 

.1103  - 

-  A 

X 

B 

C  x  D 


LHA  =  (LHA) 


15 

GHA  <  180   or   GHA  >  180 
E  =  -GHA        E  =  360-GHA 


15 


6D  =  A_ 
U         15 


15 


+ 

— 

12.0000 

XXXXXXX 

E 

W 
6D 

E 

LHA 

XXXXXXX 

TS 


(L  and  6  are  positive  if  North ,    and  negative  if  South) 


FIGURE  19.   Sunset  Calculation  Form. 
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SUNRISE 
DATE 


X 


GHA(1200) 


SIN  L 

A 

COS  L 

C 

SIN  <$ 

B 

COS  6 

D 

COS (LHA)=-. 1276  • 

-  A  x  B 

C  x  D 


LHA  =  (LHA) 


15 

GHA  <  180   or   GHA  >  180 
E  =  -GHA        E  =  36  0-GHA 


15 


6D  =  i- 
U        15 


15 


+ 

— 

12.0000 

XXXXXXX 

E 

W 
6D 

E 

LHA  XXXXXXX 

L_ 

6 


GHA(1200) 


SIN  L 

A 

COS  L 

C 

SIN  6 

B 

COS  6 

D 

COS (LHA) =- 

.1276  - 

-  A 

X 

B 

C  x  D 


LHA  =  (LHA) 


15 

GHA  <  180   or   GHA  >  180 
E  =  -GHA        E  =  360-GHA 


15 


6D  =  K 

U         15 


15 


+ 

— 

12.0000 

XXXXXXX 

E 

W 
6D 

E 

LHA  XXXXXXX 

TR 


TR 


(L  and  6  are  positive  if  North,  and  negative  if  South) 


FIGURE  20.   Sunrise  Calculation  Form. 
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The  distance  to  a  transmitting  station  should  be  greater  than 
2400  nautical  miles  to  insure  that  the  receiver  is  tracking 
the  principal  mode  of  the  transmitting  station.   The  Time 
Service  Announcement  of  the  U.  S.  Naval  Observatory,  a  sample 
of  which  is  reproduced  as  Table  IV,  is  a  necessary  document 
for  the  navigator's  use. 

It  is  also  necessary  that  the  navigator  know  the  amount 
of  frequency  offset  between  his  local  standard  and  reference 
oscillators  for  the  stations  to  be  tracked.   This  offset 
should  be  expressed  in  microseconds  and  must  be  applied  daily 
to  the  output  of  the  receivers.   As  an  exampled  of  this 
correction,  suppose  that  the  transmitter  reference  oscillator 
is  ahead  of  the  local  standard  by  1  part  in  10   .   This 
corresponds  to  a  daily  correction  of  1.16  microseconds  to  be 
added  to  the  output  data  (phase  difference)  of  the  particular 
receiver  tracking  that  transmitter.   Similar  appropriate 
corrections  will  be  applied  to  the  other  receiver  outputs. 

After  the  stations  have  been  chosen,  the  navigator  will 
check  the. antenna  field  intensity  radiation  patterns,  and 
commence  receiver  tracking.   As  the  ship  gets  underway,  the 
necessary  information  is  entered  in  the  Reference  Position 
Form.   As  time  progresses,  appropriate  entires  are  made  in  the 
Hourly  Navigation  Record.   As  stated  earlier  the   value  of 
the  oscillator  drift  was  recorded  on  the  Reference  Position 
Form.   Care  must  be  taken  to  insure  that  this  correction  and 
that  for  frequency  offset  are  properly  applied. 


66 


<*^ 

■P 

-p 

CO 

M 

P 

5                               Q)   3  * 

KD 

10 

CO    0 

X) 

• 

U          >i               -P    CO  &      • 

<J\  CO 

c 

U   WH 

N 

o  >h  4->    •  fd  fd       n 

rH  rH 

0 

3          QJ 

o 

►  -H 

C-PfJjaEH        XJoK 

rH  -H 

O   |2     CO 

CO 

^ 

0          >   OD   Ul^OAi 

U      • 

rd    CO 

^  u  fd 

• 

in  (U  u       0 

CD  O 

•H    CO 

Si 

r^ 

o 

co  in         X  o   0     •  •*  in 

•Is 

U  -H 

0    >iCu 

rH 

CO 

rH           OQJOPCOCMO 

QJ    g 

£  43 

• 

rd    g  -P        •*    a    Q) 

CU 

CU  w 

-p 

fi 

r- 

sonhfN-Hgc:^ 

•P 

CO    C 

T3    >H 

O   r- 

Cn  M   C   3        -P-hOcn 

& 

rd 

M    Q)    3 

■H   4H  -H    0    O    C   -P 

cr> 

CU 

m 

O   £   0 

d) 

c 

CO          £  £1  -P    0         QJ    C 

V£> 

CO 

Eh 

»H    O  ^ 

rH 

o 

QJ                     U    S-l  rH    0 

CT* 

T~\ 

Si 

QJ    CO  Si    d  LO          QJ  ,Q 

rH 

r^ 

i<    rH     Q) 

rd 

-p 

g   M-P  QjmNix:  fd  -P 

CM 

CO    0    C 

-P 

o 

•H    OO    >f">CO-P-P    0 

>H 

fa  ^    O 

CO 

a 

EhS^qjcmEhOcoc; 

CU 

IH 

45 

0 

O 

CO 

>i 

>1 

-P 

• 

53 

<D 

fd 

rd 

U 

a 

O 

U 

o  *t3 

O  T3 

O 

s 

H  CO 

c 

O  -H 

o  c 

CO  53 

fd 

00    U 

co   0 

<e* 

CO  O 

fi 

rH    fa 

H  S 

rH 

• 

H  H 

CU 

, 

CO 

S   Eh 
CO  < 

-P 

O  & 

-P   O 

0  Si 

■P   u 

t 

!3  H 

•H 

fd 

fd 

>i  o 

En 

3  CO 

fd 

O    CU 

O    QJ 

Pi  C* 

Pi 

S 

o 

o 

O  ro 

to 

Eh  O 

■^  Eh 

CN   Eh 

Eh  O 

qj 

H 

rH   D 

H  D 

rtj  CN 

ro 

>H    Q 

£ 

CD 

u  < 

^-~. 

Pi   . 

10 

53  K 

■O  ^ 

.— ^ 

H  U 

U 

W 

rH     Q)    ^ 

rH 

:n 

CO 

QJ 

aS 

fd  -P  — 

" — ' 

P3     • 

Qa 

c  fd 

O  Q 

3 

ss 

■H  -H    Jh 

O 

o 

CO 

g-d  cj 

O 

m 

rH       - 

— 

fa  !3 

o  fd  ^ 

o 

H 

<C  53 

2  pi  o 

** 

>  O 

Q     • 

p^ 

rH 

<  Eh 

2  CO 

•  • 

13  U 

<J 

-P 

53 

• 

u 

>. 

•  H 

W  D 

QJ 

u  * 

w  a 

S 

Uh 

c  — 

o 

o 

CO 

M 

H  g 

HH 

QJ    N 

00 

o 

.  a; 

Eh  O 

CU 

3  ffi 

• 

• 

D  & 

CO 

Pi 

&x 

r^ 

^ 

w 

fa  fa 

C 

QJ  — ' 

rH 

CM 

H 

O 

•H 

}H 

Pi 

fa 

fa 

w 

3£ 

CO 

CO 

•rH 

D 

£ 

£ 

%, 

a  z 

CD 

CT> 

a\ 

EH 

w  o 

rH 

m  • 

»  • 

53 

u 

CO 

-a 

Q) 

rH 
QJ  o 

CO 

ro 
o 

W 

X! 

C  r- 

a\ 

u 

U 

•H  V£> 

r^ 

g 

CO 

fd 

0 

o 

Cn 

c 

!3 

0  S 

a 

•H 

0 
•H 

5h  -• 

•>  Ni  ro 

rd 

3 

g 

-P 
fd 

Q)   00 
rH  ro 

0  rH  ro 
42    fd  o 

w 

rH 

u 

-P  o 

rH      C    O 

u 

iH 

o 

3  ^T 

rd    rd  cr> 

H 

O 

HH 

U  ^r 

fflUo 

> 

MH 

Pi 

w 

CU 

a 

w 

4h 

EH 

o 

•H 

■p 

fd 

2 

<3 
PQ 

H 

• 

•p 

B 

2 

EH 

rH 

CO 

67 


to 
O 

rH  -H 

rd  co 
•H  CO 
U  -H 
0  g 
ft  CO 
CO  £ 
rd 
H 

EH 


rd  z 
c  fd 

•H  -H    $H 

o  d  :-j 

2^0 


Pk 


CD    N 

CD  *-- 

U 

En 


(1) 


1 

,c 

rH 

C 

u 

-Q 

to 

•H 

rd  T) 

fd 

3 

e 

<D 

Q) 

p 

0 

X 

co 

3 

0) 

<D 

O 

c 

> 

in 

0 

0) 

•H 

•H 

0 

r-\ 

to 

4-) 

U-l 

IM 

fd 

c 

CD 

c  x: 

0 

P 

X5 

0  Oi 

u 

ft 

CD 

CO 

u 

• 

w 

u 

CD 

2 

i>i 

co 

X 

p 

c 

<D 

Em 

CD 

3  x:  x 

Eh 

>i 

D 

rd 

'd  ^ 

o 

co 

p 

o 

u 

a 

LD 

3 

0 

H  ,£ 

g 

En 

O 

Xi 

P  T 

u 

£ 

td 

O 

0 

a) 

o 

U 

o 

CD 

MH 

H 

CO 

O 

-p 
o 

P    N 

O  LT) 

>£>  »£> 

00  00 
rH   H 

c  c 

o  o 


o 

CM 


o 


CO 
rH 


^ 

to 

rrj 

r--_ 

t^» 

»* 

^ 

k 

£ 

0) 

j?j 

o 

CD 

CM 

•-  • 

H 

» . 

Lfl 

U  (N 

in 

rH 

0 

g 

0 

H 

•H 

CO 

CM 

h) 

^r 

rH 

2 


En  Pl, 


>i 

rd 
■d 
c 
o 
a 


Eh 
D 

o 

o 

rH 


O 


C 

o 

p 
o 
c  * 

p    N 

<-'      H-* 

X!  ^ 

O  LD 

ro  ro 

CN  CM 


o 

CO    >,,£ 


tH  fd 

>i  C 

fd  o 

CO 

CD  T3 

P  rH 

Eh  ro 


O 


ro 

CM 


•H 

fd 

Z 

fd 

ffi 

K 

•H 

£ 

0)  53 

ro 

rH  m 

—  • 

rd  -• 

en 

p  m 

o 

rH   CM 

o 

fd  o 

CO 

3  -H 

in 

(^    CM    rH 

a 

53 


0  P 
H    3 

•H 

e 


w 


fd  ,c 

C  P 
tr.  o 
•H  ^o 
in 

0 
0)  P 

e 

•H  in 

Eh  m 


U 


H  P 
0 


CO 

0  • 

3  CD 

C  rH 

H  & 

rd 

r:  p 

o  w 


u 

0) 
CD    CO 


o  s: 
o  u 
ro  rd 
rH    rj 


in 

CO  ( 


o 


CM 


CM 
CM 

a  vr> 

co    ■> 
•H  a 

rH    rH 

0   -• 

rd  m 
C  o 
C  oo 
rtj  ro 


f73 
2 


Eh 

rd 

O  ^ 

o  CO 

en  m 

rH  C 

0  CD 

P  £ 


CD 
CO 
rd 
XX 
0* 


U 

C 
rd 

CO 

Ed 


in 

o  ro 
ro     . 

•  CM 
CM  CM 
CM 
C 
0 

P 

o 
c 

p 
3 


^r 


^3  X 


Eh  — 
p  ro 

o 
o    >i 

ro   rd 

O    ^ 

c 


0 
P 


o  x: 
o  u 
o  rd 
o   QJ 


o 
o 
o 


rH 


o 
ro 


CM 
CM 


tt) 

ft 

rd 

U 

P 

10 

rd 

V 

w 

CD 

•H 

w 

CO 

&   r 

o 

rd 

-  • 

en 

XI 

H 

o> 

o 

p 

p 

■^ 

0 

H 

co 

0 

^r 

0 

p 

rH 

r-i 

53 

rtj 

CM 

rH 

u 


p 
c. 

I 

0) 

u 

c 

o 


0) 
U 
•H 
> 
U 
0) 
CO 

I 

•H 

Eh 


> 

H 

W 
in* 

9 


68 


The  navigator  will  calculate  the  time  of  sunset  for  each 

of  the  transmitting  stations  using  a  Sunset  Calculation  Form. 

Using  the  ship's  speed  and  course,  he  will  estimate  his 

position  at  ionospheric  sunset,  and  calculate  this  time  of 

sunset.   These  calculation  will  require  the  use  of  the  Nautical 

Almanac  to  determine  the  declination  and  the  1200  Greenwich 

Hour  Angle  for  the  sun.   Now  using  a  modified  version  of 

Wait's  equation  for  the  amount  of  diurnal  shift  (the  equation 

now  allows  the  use  of  distance  in  nautical  miles)  we  arrive 

at  the  value  of  At. 

At  =  D  C 
v 

The  value  of  D  is  readily  accessible  from  the  tables  of  posi- 
tion.  Again  using  predicted  course  and  speed,  the  various 
times  of  sunrise  may  be  computed.   C   is  obtained  from  Figures 
4  and/or  6,  depending  on  the  propagation  path.   Sufficient 
information  is  now  at  hand  for  the  diurnal  shift  prediction 
trapezoid.   The  trapezoid  must  be  modified  to  indicate  a  slope 
due  to  the  motion  of  the  ship  with  respect  to  the  transmitting 
station.   The  amount  of  correction  for  each  hour  is  N  times 
6.18  microseconds,  where  N  is  the  ship's  speed  component  in 
knots  toward  (minus  correction)  or  away  from  (plus  correction) 
the  transmitter.   With  the  construction  of  the  modified 
trapezoid  completed,  the  hourly  corrections  are  entered  on 
the  bottom  of  the  Hourly  Navigation  Record  for  use  at  the 
required  time. 

As  an  illustration  of  the  accuracy  of  this  method  for  the 
prediction  of  At,  Table  V  lists  the  following  results  for  two 

stations,  NPM  and  NAA. 
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TABLE  V 
TRAPEZOID  PREDICTION  RESULTS 

NPM 

frequency  -  23.4  kHz 

propagation  path  -  seawater 

predicted  At  -  41.2  microseconds 

mean  error  -  2.5  microseconds  (810  yards) 

NAA 

frequency  -  17.8  kHz 

propagation  path  -  land 

predicted  At  -  30.3  microseconds 

mean  error  -  2.1  microseconds  (680  yards) 
In  the  above  two  examples,  mean  error  is  equal  to  the  average 
of  the  absolute  values  of  the  differences  between  the  predicted 
trapezoid  and  the  recorded  diurnal  shift.   Samples  of  the 
differences  were  taken  every  fifteen  minutes  during  the  period 
that  diurnal  shift  was  present. 

Occasionally,  it  is  necessary  for  transmitting  stations 
to  make  offsets  to  their  master  oscillators  to  stay  in  tolerance 
with  the  national  standard.   These  offsets  are  scheduled  well 
in  advance  of  their  accomplishment.   The  navigator  must  make 
a  similar  correction  to  his  readings  at  the  same  time  as  the 
actual  offset  is  applied  to  maintain  a  tracking  navigation 
system.   This  offset  is  applied  only  if  navigation  is  in 
process.   This  correction  can  be  ignored  when  starting  from 
a  new  reference  position,  since  in  this  case,  the  offset  is 
then  included  in  the  A  value. 
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APPENDIX  A 
STATIONS 

NPM  NAA 

L  =  21-25. 5N  L  =  44-38. 9N 

X  =158-09. 3W  X  =  67-16. 9W 

158.1550W  67.2817W 

A  =  .3652  A  =  .7028 

C  =  .9308  C  =  .7114 

0D  =  10.5437W  eD  =  4.4854W 


NBA 

NSS 

NWC 

L  = 

09-03. 3N 

L 

= 

38-59. IN 

L 

=  21-49. OS 

X  = 

79-38. 9W 
79.6483W 

X 

= 

76-27. 2W 
76.4533W 

X 

=114-09. 8E 
114.1633E 

A  = 

.1574 

A 

= 

.6291 

A 

=-.3716 

C  = 

.9876 

C 

= 

.7773 

C 

=  .9284 

6D  = 

5.3099W 

SD 

= 

5.0969W 

9D 

=  7.6109E 

GBR 

NLK 

WWVL 

L  =  52-22N 

L 

=  48-12N 

L  =  40-40N 

X  =   1-11W 

X 

=121-55W 

X  =105-02W 

1.1833W 

121.9167W 

105.0333W 

A  =  .7919 

A 

=  .7455 

A  =  .6517 

C  =  .6106 

C 

=  .6665 

C  =  .7585 

6D  =  .0789W 

6D 

=  8.1278W 

0  -=  7.0022W 
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Time  or  X 

in  Decimals 

.0000 

.0083 

.0167 

.0250 

.0333 

.0416 

.0500 

.0583 

.0667 

.0750 

.0833 

.0916 

.1000 

.1083 

.1167 

.1250  ' 

.1333 

.1416 

.1500 

.1583 

.1667 

.1750 

.1833 

.1916 

.2000 

.2083 

.2167 

.2500 

.2333 

.2416 

.2500 

.2583 

.2667 

.2750 

Time  or  A 
in  Minutes 

00 


01 
02 
03 
04 
05 
06 
07 
08 
09 
10 
11 
12 
13 
14 
15 
16 


Time  or  X    Time  or  X 
in  Decimals   in  Minutes 


.0017 

.1 

.0033 

.2 

.0050 

.3 

.0067 

.4 

.0083 

.5 

.0100 

.6 

.0117 

.7 

.0133 

.8 

.0150 

.9 
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Time  or  X 

Time  or  A 

in  Decimals 

in  Minutes 

.2833 

17 

.2916 

.3000 

18 

.3083 

.3167 

19 

.3250 

.3333 

20 

.3416 

.3500 

21 

.3583 

.3667 

22 

.3750 

.3833 

23 

.3916 

.4000 

24 

.4083 

.4167 

25 

.4250 

.4333 

26 

.4416 

.4500 

27 

.4583 

.4667 

28 

.4750 

.4833 

29 

.4916 

.5000 

30 

.5083 

.5167 

31 

.5250 

.5333 

32 

.5416 

.5500 

33 

.5583 
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Time  or  X 

Time  or  X 

in  Decimals 

in  Minutes 

.5667 

34 

.5750 

.5833 

35 

.5916 

.6000 

36 

.6083 

.6167 

37 

.6250 

.6333 

38 

.6416 

.6500 

39 

.6583 

.6667 

40 

.6750 

.6833 

41 

.6916 

.7000 

42 

.7083 

.7167 

43 

.7250 

.7333 

44 

.7416 

.7500 

45 

.7583 

.7667 

46 

.7750 

.7833 

47 

.7916 

.8000 

48 

.8083 

.8167 

49 

.8250 

.8333 

50 

.8416 
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Time  or  X 

Time  or  X 

in  Decimals 

in  Minutes 

.8500 

51 

.8583 

.8667 

52 

.8750 

.8833 

53 

.8916 

.9000 

54 

.9083 

.9167 

55 

.9250 

.9333 

56 

.9416 

.9500 

57 

.9583 

.9667 

58 

.9750  ' 

.9833 

59 

.9916 

1.0000 

60 
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NPM         23.4  KHZ  21-25N    158-09W 

LATITUDE  LONGITUDE  WEST 

NORTH      114-00  114-15  114-30  114-45 

00-00      2885.37  2872.30  2859.26  2846.22 

17831.57  17750.84  17670.20  17589.65 

299.4/029.4  299.5/029.5  299.6/029.6  299.7/029.7 

OC-15      2878.03  2864.94  2851.87  -2838.81 

17786.22  17705.34  17624.54  17543.83 

299.2/029.2  299.3/029.3  299.4/029.4  299.5/029.5 

00-30      2870.74  2857.62  2844.52  2831.44 

17741.15  17660.11  17579.16  17498.30 

298.0/029.0  299.1/029.1  299.2/029.2  299.3/029.3 

00-45      2863.49  2850.35  2837.23  2824.12 

17696.36  17615.17  17534.07  17453.04 

298.8/028.8  298.9/028.9  299.0/029.0  299.1/029.1 

01-00      2856.29  2843.13  2829.98  2816.84 

17651.86  17570.52  17489.26  17408.08 
298.6/028.6  298.7/028.7  298.8/028.8  298.9/028.9 

01-15      2849.13  2835.95  2822.77  2809.61 

176C7.63  17526.14  17444.74  17363.41 

298.4/028.4  298.5/028.5  298.6/028.6  298.7/028.7 

01-30      2842.02  2828.81  2815.62  2802.43 

17563.70  17432.06  17400.50  17319.03 

298.2/028.2  298.3/028.3  298.4/028.4  298.5/028.5 

01-45      2834.96  2821.73  2808.51  2795.30 

17520.05  17438.27  17356.57  17274.94 

298.0/028.0  298.1/028.1  298.2/028.2  298.3/028.3 

02-00      2827.94  2814.69  2801.44  2788.21 

17476.70  17394.77  17312.92  17231.15 

297.8/027.8  297.9/027.9  298.0/028.0  298.1/028.1 

02-15      2820.98  2807.70  2794.43  2781.17 

17433.64  17351.57  17269.58  17187.66 

297.6/027.6  297.7/027.7  297.8/027.8  297.9/027.9 

02-30      2814.06  2800.75  2787.46  2774.19 

17390.87  17308.67  17226.53  17144.47 
297.4/027.4  297.5/027.5  297.6/027.6  297.7/027.7 

02-45      2807.19  2793.86  2780.55  2767.25 

17348.41  17266.06  17183.78  17101.58 

297.2/027.2  297.3/027.3  297.4/027.4  297.5/027.5 
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10 


11 

13 


15 
16 


25 

28 

26 
27 


OMPUTES  TABLES  OF  POSITIONS 

36,8),  0BM(36,S) 

TS,  ALNGS 

) 

T,  ALNG 


THIS  PROGRAM  C 

DIMENSICN  CPM( 

READ  (5,3)  ALA 

FORMAT  (2F10.5 

READ  (5,3)  ALA 

AMALT  =  40. 

ALNG1  =  ALNG 

IAB  =  34 

IF  (ALAT.GE. AMALT)   GO  TO  650 

II  =  1 


DO  500 

DO  400 

RALAT  = 

RALATS 

IF 

IF 


I  =  It 

J  =  1, 

ALAT/5 

ALATS 

(ALNGS.GE.O 

(ALNG.LT.O. 


DLO  =  ALNGS  - 

ACLO  =  ABS(DLO 

IF  (ADL0.GE.18 

GO  TO  13 

IF  (ALNG.GE.O. 

DLO  =  ALNGS  - 

IF  (DLO. GE. 180 

GO  TO  13 

DLO  -    ALNGS  - 

GG  TO  13 

ADLO  =  ABS(DLO 

CONTINUE 

HADLO  =  C.5-AD 

THE  =  SIN(HADL 

DL  =  ALATS  -  A 

ADL  =  ABS(DL) 

HADL  =  0.5*ADL 

HD  =  SIN(HADL) 

CD  =  1.  -  2.-*H 

IF  (CD.GT.O.O) 

B  =  -CD 

ARGS  =  B/SQRT( 

RD  =  1.5708  + 

GO  TO  16 

ARGC  =  CD/SQRT 

RD  =  1.5708  - 

CONTINUE 

D  =  RD  *  57.29 

DNM  =  60.  *  D 

DMS  =  6.18  *  D 

CCL1  =  90.  -  A 

AD  =  D  -  CCL1 

DCL1  =  ABS(AD) 

IF  (ALAT.GE. 0. 

IF  (ALATS. GE.O 

GO  TO  26 

IF  (ALATS. GE.O 

CCL2  =  90.  +  A 

GO  TO  27 

C0L2  =  90. 

CONTINUE 

HC0L2  =  0.5' 

HCCL1  =  C.5; 


12 

8 
7.29578 
/57. 29578 
.0)   GO  TO  10 
0)   GO  TO  11 
ALNG 
) 
0. )   DLO  =  360 

0)   GO  TG  11 

ALNG 

. )   DLO  =  -360 


+  DLO 


+  DLO 


ALNG 

) 

LO/57. 29578 

0 ) *S  IN (HADLO )*COS( RALAT )*COS( RALATS) 

LAT 

/57. 29578 
•SIN(HADL)     +    THE 
D 

GO  TO  15 

1.  -  B*^2) 
ATAN( ARGS) 

(  1.  -  CD-s^2) 
ATAN(ARGC) 

578 

NM 
BS(ALAT) 


0)   GO  TO  25 
.0)   GO  TO  28 

.0)   GO  TO  26 
BS(ALATS) 


-  ABS(ALATS) 


CA  = 
CB  = 
CC  =  1. 
HCCC  = 
CCCC  = 
ARGCC  = 


CO 
DC 


SIN(HCGL2 


/COS(RA 
/SIN(RD 
CA  *  CB 
1.  -  2. 

cccc/s 


30 


RCCC  =  1.5708 

CCC  =  RCCC  *  5 

IF  (ALAT.GE.O. 

IF  (DLO.GE.C.O 

GO  TO  36 

IF  (DLO. GE. 0.0)   GO  TO  40 

GC  TO  37 


L2/57. 29578 
Ll/57. 29578 

)vSIN(HC0L2)  -  SIN(HDCL1)*SIN(HDCL1) 
LAT) 
) 

#  CC 
-:  HCC C 

QRTd.    -    CCCC*&2) 
-    ATAN(ARGCC) 
7.29578 
0)       GO    TO    30 
)       GO    TO    35 
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35  CCC    =    CCC    -    180 
GO    TO    40 

36  CCC    =    CCC    +    180. 
GO    TO    40 

37  CCC  =  360.  -  CCC 
40  CONTINUE 

BCC  =  CCC  +  90. 

IF  (BCC. GT. 360. )   BCC  =  BCC  -  360. 

IJ  =  I  I  +  1 

IK  =  II  +  2 

OPM( II  ,  J)  =  DNM 

OPM(IJ,J)  =  CMS 

OPM(  IK,  J)  =  CCC 

OBM< IK, J)  =  BCC 

ALNG  =  ALNG  -  0.25 
4C0  CCNTINUE 

ALNG  =  ALNG1 

I  I  =  IK  +  1 

ALAT  =  ALAT  +  0.25 
500  CCNTINUE 

WRITE  (6,100) 

100  FORMAT  (1H1,"NWC",7X,"22.3  KHZ" , 2X , "2 1-49 .OS" , 4X , 
1"114-09.8E") 

WRITE  (6,101) 

101  FORMAT  ("LATITUDE", 47X, "LONGITUDE  WEST",47X, 
1"LATITUDE") 

WRITE  (6,102) 

102  FORMAT  ( 2X, "NORTH" , 5X , "124-00" , 7X, " 124-1 5 "  ,7X , 
1"124-30",7X,"124-45",7X,"12  5-00",7X,"12  5-15",7X, 
1 "125-30", 7 X,"l 2  5-45", 9X, "NORTH") 

L  =  1 

ICD  =  00 

DO  600   K  =  1,  36 

GO  TO  (601,602, 603)  L 
6C1  WRITE  (6,605)  I AB, I  CD, OPM ( K , 1 ) ,0PM ( K , 2 ) ,0PM ( K , 3 ) , 

10PM(K,4),0PM(K,5),0PM(K,6),0PM(K,7),0PM(K,8),IAB,ICD 
6C5  FORMAT  ( 1H0 , IX , I  2 , "-",  12 , 4X, F8 . 2, 5X, F 8. 2 , 5X , F8. 2 , 5X, 
1F8.2,5X,F8.2,5X,F8.2,5X,F8.2,5X,F8.2,8X, 12,"-", 12) 

ICD  =  ICD  +  15 

IF  (ICD.EQ.60)   GO  TO  615 

GO  TO  610 
615  ICD  =  00 

IAB  =  IAB  +  1 

GO  TO  610 

602  WRITE  (6,606)  OPM( K , 1 )  , OPM ( K , 2  )  ,OPN ( K ,3 )  , OPM ( K,4 )  , 
10PM(K,5),0PM(K,6),0PM(K,7),0PM(K,8) 

606  FORMAT  ( 1 IX , 8( F 8 .2 , 5X) ) 
GC  TO  610 

603  WRITE  (6,607)  OPM( K , 1 ) , OBM ( K , 1 ) , OPM ( K , 2 ) , OBM ( K , 2 ) , 
10PM (K, 3) ,0BM(K,3) ,0PM(K,4) ,0BM(K,4) ,0PM(K,5) ,0BM(K,5) , 
10PM(K,6) ,03M(K,6),0PM(K,7)  , OBM ( K , 7 ) , OPM ( K  ,  8 ) ,CBM(K,8) 

6C7  FORMAT  ( <3X  ,  8  (  F  5.  1 ,  "  /"  ,  F5  . 1 ,  2X  )  ) 

GC  TO  610 
610  L  =  L  +  1 

IF  (L.GT.3)   L  =  L/3 
600  CONTINUE 

GO  TO  7 

650  WRITE  (6,651) 

651  FORMAT  (1H1) 
END 
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